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CHAPTER  I 
INTRODUCTION 


1 . 1  Background 

Mesospheric  and  stratospheric  water  vapor  have 
historically  been  elusive  quantities  to  measure,  and  yet 
such  measurements  are  of  prime  importance  in  understanding 
atmospheric  processes.  Water  vapor  is  a  key  constituent 
in  many  atmospheric  chemical  reactions,  and  is  of  obvious 
importance  to  ice  nucleation  processes  with  their  conse¬ 
quent  release  of  latent  heat.  A  variety  of  atmospheric 
models  have  been  developed  with  widely  varying  water  vapor 
mixing  ratios.  Measurements,  moreover,  have  been  plagued 
with  uncertainties  since  the  earliest  attempts  at  measure¬ 
ment  in  the  nineteen  forties . 

The  first  in-situ  measurements  were  suspected  of 
contamination  by  the  measuring  instruments  and  balloons. 
During  the  late  nineteen  fifties,  methods  of  remotely 
measuring  water  vapor  were  attempted.  Remote  measurements 
were  free  of  contamination,  however,  the  increase  difficulty 
encountered  in  conducting  remote  measurements  lowered  their 
credibility.  Current  remote  measurements  vary  as  much  as 
75  ppm  whereas  current  models  allow  only  about  10  ppm 
variation.  Thus  the  measurement  of  high  altitude  water 
vapor  is  still  in  the  research  phase  and  should  improve 
greatly  with  advances  in  theory  and  technology. 


the  Problem 

Remote  water  vapor  measurements  take  advantage  of 
the  H.,0  molecule's  natural  resonant  properties  in  the 
microwave  portion  of  the  frequency  spectrum.  Figure  1 
(Longbothum,  1976)  shows  three  resonances  of  water  vapor 
and  how  atmospheric  oxygen  raises  the  apparent  background 
radiation  about  each  line.  Figure  2  (Longbothum,  1976) 
shows  the  apparent  pressure  broadening  occurring  in  the 
lower  altitudes  for  the  22  and  183  GHz  lines.  The  line 
shape  reveals  not  only  an  approximate  water  vapor  mixing 
ratio,  but  may  provide  mixing  ratio  versus  altitude  data 
if  suitable  data  analysis  is  employed. 

Two  types  of  measurements  are  considered  here, 
emission  and  absorption.  Resonant  emissions  are  detected 
by  observing  the  atmosphere  against  a  background,  as  free 
as  possible  of  other  sources,  by  looking  away  from  galacti 
sources.  Resonant  absorption  is  detected  by  observing  a 
hot  source  beyond  the  atmosphere,  most  generally  the  sun. 
The  intervening  water  vapor  will  then  absorb  radiation 
according  to  its  mixing  ratio  and  the  relative  atmospheric 
pressure . 

It  was  decided  to  detect  the  water  vapor  resonance 
centered  at  22.23507985  GHz,  since  this  is  the  only  line 
which  may  be  studied  from  the  ground.  Although  the  peak 
amplitude  at  22  GHz  is  relatively  small,  according  to 
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Figure  1,  the  overall  atmospheric  attenuation  is  also  small. 
At  183  GHz  and  above,  the  atmospheric  absorption  is  so 
great  that  the  higher  altitude  contribution  would  be 
unusable  to  a  ground  based  detector.  Other  attractive 
features  of  the  22  GHz  line  are  its  relative  remoteness 
from  other  lines  and  its  relatively  low  frequency.  The 
remoteness  will  cause  least  bias  from  other  resonances, 
and  the  lower  frequency  employs  components  which  are  less 
expensive  and  more  readily  available. 

An  instrument  capable  of  detecting  the  22  GHz  line 
was  proposed  by  Longbothum  (1976)  and  has  since  been 
improved  upon.  Longbothum  describes  a  microwave  radio¬ 
meter  (Figure  3)  consisting  of  a  super  heterodyne  receiver 
feeding  a  large  filter  bank.  Each  filter  then  feeds  a 
detector,  integrator,  and  analog  to  digital  converter. 

A  computer  would  then  periodically  sample  the  A  to  D 
outputs  and  continue  to  compile  and  finally  reduce  the 
data . 

In  recent  times,  autocorrelators  have  been  exten¬ 
sively  employed  as  spectral  analyzers,  especially  for  radio 
astronomy  applications.  Autocorrelators  have  inherent 
advantages  over  filter  banks  as  will  be  discussed  herein. 
Thus ,  an  improvement  over  Longbothum' s  scheme  is  to 
replace  the  filter  bank  and  associated  circuitry  with  an 
autocorrelator.  The  implications  and  analysis  of  this 
change  will  be  described  in  the  following  chapters. 


I  ■  . . .  i.. .  - 
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CHAPTER  II 
DESIGN  CRITERIA 


2 . I  Basic  System  Requirements 

A  system  is  needed  that  will  detect  the  spectral 
content  of  the  water  vapor  line  centered  at  22.23507985  GHz. 
The  desired  altitude  of  the  line  will  be  from  approximately 
40  km  up . 

Figure  2  shows  that  the  line  halfwidth  corresponding 
to  40  km  is  10  MHz.  The  necessary  system  bandwidth  then 
becomes  20  MHz.  This  bandpass  effectively  reduces  the 
lower  altitude  information  to  an  approximate  DC  bias  on 
the  high  altitude  curve.  This  permits  a  high  altitude 
measurement  from  a  ground  based  observation. 

Although  this  report  will  not  be  concerned  with  RF 
hardware,  a  brief  description  of  overall  system  sensitiv¬ 
ity  and  calibration  will  be  required.  Figures  1  and  5 
are  a  theoretical  representation  of  the  water  vapor  emission 
and  absorption  spectra,  respectively.  These  have  been  de¬ 
rived  from  a  computer  simulated,  ground  based  observation 
of  the  atmosphere  (Longbothum,  1977) .  As  was  mentioned 
previously,  water  vapor  measurements  vary  greatly  and  thus 
the  atmospheric  models  vary  also.  These  simulations  are 
based  on  only  one  such  model  and  should  therefore  be  con¬ 
sidered  only  as  guidelines  for  a  svstem  design.  This 
model  and  others  are  more  fully  described  in  Longbothum' s 


■  ^  ••  -  • 


report . 


on  spectra  (Longbothum,  1977) 


These  simulations,  however,  do 


em  sensitivitie 


emission  experiment,  temperatures 


of  173  K  must  be  received  with  a  system  sensitivity  to  sig 


nals  0.5  K  or  less.  This  sensitivity  will  be  defined  as 


the  minimum  detectable  signal,  IT  .  .  Although  IT  .  is 

°  mm  °  mm 

given  here  as  the  signal  amplitude,  in  practice  a  much 
smaller  &T  ^  will  be  required  to  provide  any  accuracy  in 
line  shape. 


'he  requirements  of  the  absorption  experiment  are 


equiring  antenna  temperatures  of  4300  K  with 


a  minimum  detectable 


To  improve  measurement  accuracy,  that  is,  to  reduce 


measurement  variance,  the  minimum  detectable  signal  should 
be  made  as  small  as  possible.  These  extreme  sensitivities 
can  be  accomplished  by  accumulating  the  desired  signal  in 
time.  The  following  relationship  defines  the  radiometer's 


minimum  detectable  signal  when  viewing  a 


(Adapted  from  Tiuri,  1964,  and  O' Kean  and  Lombardo 


aTon 


min 


li 


’.'.'here : 

it  =  Coefficient  describing  type  of  receiver 
BW  ■  High  frequency  bandwidth 
Top  =  Svstcm  operating  temperature 
tt  =  Integration  time 

The  AT  .  here  is  equivalent  to  the  standard  deviation  of 
nun  ^ 

the  received  power. 

When  viewing  a  thermal  source  over  time,  amplifier 
gains  mav  vary  somewhat,  making  a  calibration  step  essen¬ 
tial.  The  Dicke  radiometer  (Dicke,  1946)  mechanically 
switches  the  radiometer  to  a  standard  noise  source  30 
times  per  second.  This  rapid  transition  gives  a  reference 
signal  often  enough  to  reduce  gain  variations.  An  alter¬ 
native  to  Dicke  switching  is  frequency  switching.  In  this 
method  the  local  oscillator  frequency  is  switched  such  that 
the  radiometer  is  looking  off  line  to  a  flat  portion  of 
the  spectrum.  Then,  by  accumulating  the  difference  between 
the  on  and  off  frequency  measurements  a  relative  spectrum 
may  be  determined,  independent  of  receiver  or  signal  DC 
characteristics.  Then,  to  achieve  a  true  record  of  receiver 
gain  variations,  a  waveguide  switch  is  used  to  switch  from 
the  antenna  to  a  reference  noise  source.  These  two  refer¬ 
ences  and  the  signal  measurement  must  then  be  taken  before 
the  receiver  gain  varies.  Thus,  the  switching  time  must 
be  short  but  also  long  enough  to  make  the  receiver  local 
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oscillator  lock  up  time  and  the  waveguide  switching  time 
insignificant.  A  reasonable  measurement  period  would  be 
one  second  in  each  mode.  The  spectral  analyzer  must  then 
cope  with  the  various  calibrating  signals  over  the  time 
necessary  to  achieve  the  desired  sensitivity. 

The  accuracy  of  the  water  vapor  mixing  ratios  is 
dependent  on  the  measured  spectral  resolution.  That  is, 
each  spectral  point  is  a  function  of  frequency  and  altitude, 
and  can  be  used  to  determine  the  water  vapor  mixing  ratio  at 
that  particular  altitude.  This  is  accomplished  by  solving 
the  radioactive  transfer  equation  developed  by  Chandrasekhar 
(1960).  Determining  the  actual  water  vapor  mixing  ratios 
is  beyond  the  scope  of  this  paper  but  is  given  some  con¬ 
sideration  by  Waters  (1976)  and  Longbothum  (1976) .  Spec¬ 
tral  resolution,  then,  greatly  affects  the  resultant  water 
vapor  mixing  ratios  and  will  eventually  require  much  atten¬ 
tion.  The  optimum  number  of  points  required,  the  centers 
and  bandwidths  are  also  beyond  the  scope  of  this  paper. 

A  sufficient  criterion  for  the  purposes  here  is  simply  to 
provide  as  accurate  a  line  representation  as  possible. 
Included  in  that  representation  is  a  requirement  of  having 
one  spectral  point  located  at  the  exact  center  frequency 
of  the  water  vapor  line. 


In  conclusion,  the  desired  spectral  analyzer  is 
required  to  analyze  a  band  20  MHz  wide.  To  prevent 
aliasing,  this  allows  a  center  frequency  no  lower  than 
10  MHz.  The  receiver  minimum  detectable  signal  must  be 
less  than  0.5  K,  and  the  spectral  analyzer  is  required  to 
remain  stable  for  the  amount  of  time  necessary  to  achieve 
the  desired  sensitivity.  The  analyzer  then  is  required 
to  produce  a  calibrated  power  spectrum  over  the  20  MHz 
band  with  emphasis  at  the  center  frequency.  The  basic 
radiometer  now  appears  as  shown  in  Figure  6. 


CHAPTER  III 


METHODS  OF  SPECTRAL  ANALYSIS 

3.1  Fi Iter  Method 

Power  spectra,  basically,  can  be  effee'  ivelv 
measured  using  two  different  techniques.  (Another  method, 
the  swept  frequency  or  "spectrum  analyser"  technique, 
sacrifices  too  much  of  the  available  information  to  be 
seriously  considered.)  The  first  technique,  as  was 
suggested  by  Longbothum,  involves  dividing,  the  signal 
through  a  bank  of  bandpass  filters.  The  filter  outputs 
are  then  squared  and  averaged  in  time  to  provide  the 
power  spectrum.  This  method  can  be  accomplished  using 
either  digital  or  analog  circuitry,  or  combinations  of 
the  two.  The  basic  digital  and  analog  circuits  are 
shown  in  Figures  7  and  8,  respectively. 

The  analog  method  is  straightforward,  however,  in 
order  to  realize  the  scope  of  digital  requirements,  the 
digital  filter  deserves  somewhat  further  explanation. 

The  most  general  method  of  implementing  digital  filters 
is  the  recursive  realization.  That  is,  a  filter  whose 
transfer  function  includes  both  polos  and  zeros,  as 
described  here  in  terms  of  the  z-transform  (Rabiner  and 
Gold,  1975) . 


Figure  8.  Digital  filter  method  of  spectral  analysi 
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The  difference  equation  is  then  derived  using  the  inverse 
2- transform: 

N  N 

v(n)  “  I  a.x(n-i)  -  Z  b  .  v  (n-  i") 
i-0  1  i- 1  L' 

A  realisation  of  this  filter  is  shown  in  Figure  9.  The 
number  of  delays  is  analogous  to  the  order  of  an  analog 
filter,  predominantly  controlling  the  transition  band 
phase  and  the  passband  and  stopband  ripple.  Cutoff  fre¬ 
quencies  are  sec  by  the  proper  combination  of  the  sampling 
rate  and  multiplying  constants. 

A  special  case  of  the  recursive  filter  leads  to  a 
second  class  of  digital  filters,  namely  the  nonrecursive 
filters.  These  filters  are  characterised  by  having  poles 
at  z^O  only,  that  is,  the  denominator  of  Equation  (2''  is 
equal  to  1.  This  leads  to  a  difference  equation  as  shown 
in  Equation  (4),  and  a  realization  as  show,  in  Figure  10. 

N 

y(n)  =  v  a.x(n-i)  (4) 

i=0  1 

Both  of  these  filters  are  discussed  extensively  in 
a  variety  of  textbooks  (Rabiner  and  Gold,  1°75,  Oppenhein 
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and  Schafer,  1975).  The  intent  here  is  only  to  attain  a 
feel  for  the  complexity  of  the  circuitry.  An  important 
consideration  is  the  fact  that  for  a  g,iven  filter  transfer 
function,  the  likelihood  of  obtaining  integer  values  of 
a^  and  is  negligible.  This  means  that  either  a  complex 
multiplier  is  required  or  the  intended  filter  characteris¬ 
tics  must  be  sacrificed. 

3 . 2  Autocorrelation  Method 

A  second  technique  for  measuring  power  spectra  is 
autocorrelation.  The  basic  autocorrelator  performs  the 
autocorrelation  function,  R(iAr) ,  on  the  incoming  signal, 
x(t)  . 

1  fT 

R(iA-r)  =  sr  x(t)x(t+iAr)dt ;  i=0 , 1 , 2 . N-l  (5) 

1  o 


Where : 

T  =  Total  correlation  time 
Ar  =  Unit  delay 
N  =  Correlation  points 
or  digitally; 

lk-l 

R(iAt)  *  ;r  I  x(kAt)x(kAt+iAT) ; 
K  k=0 

i-0. 1,2 . N-l 


(6) 


K  =  Number  of  products 
At  =>  Sampling  interval 

The  result  is  then  multiplied  by  a  weighting  function, 
w(y) ,  and  the  Discrete  Fourier  Transform  is  performed 
on  the  products  (Weinreb,  1963). 

N-l 

P(i<5f)  =  2AtR(o)w(o)+4At  Z  R(nAT)w(nAr)  cos 

n=l 

(2TTi5fnt) 


Where : 

<Sf  =  Frequency  difference  between  filter  centers. 

The  weighting  function  is  required  to  make  the 
result  here  equal  to  the  result  of  a  discrete  filter  out¬ 
put  (Schwartz  and  Shaw,  1975,  Weinreb,  1963).  The  analog 
form  of  the  filter  output  is  shown  in  Figure  7  and  repeated 
in  Equation  (8),  while  its  equivalent  digital  form  is  shown 
in  Figure  8  and  Equation  (9)  . 


f  T 

P(6f)  =  ^  j  x2(t)dt 

'  O 

i  K-l  0 

P(6f)  =  i  Z  (kAt) 
K  k=0 


(3) 

(9) 


Weinreb  shows  an  equivalence  between  Equations  (5)  and  (7) , 
providing  the  weighting  function  is  any  function  exhibiting 
the  following  conditions: 
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(a)  w(o)  -  1  =  W(f)df  (10) 

4  —  »»' 

Inhere  W(f)  is  the  fourier  transform  of  w(o)  ; 

(b)  w(nA')  ■  w(-nAi)  (11) 

(c)  w(nAx)  »  0  for  'n  N  (121 

Inhere  N  is  the  number  of  delays  used  in  the 
autoeorrelat ion  function. 

A  number  of  weighting  functions  are  in  common  use  (Weinreb, 
1963,  Hamming,  1977,  Schwartz  and  Shaw,  1975)  ,  however, 
according  to  Blackman  and  Tukev  (1958),  the  choice  of 
window  does  not  considerably  influence  the  resulting  spec¬ 
trum  . 

A  disadvantage  in  using  autocorrel ator s  over  the 
filter  method  of  spectral  analysis  is  in  the  spacing  be¬ 
tween  frequency  points.  Weinreb  (19e3)  shows  that  evenly 
spaced  delavs  result  in  evenly  spaced  spectral  points. 

That  spacing,  3f,  corresponds  to  1/N.\r,  where  N  is  the 
number  of  autocorrelation  delays  and  At  is  the  unit  delav 
time.  The  half  power  bandwidth,  Af,  of  each  spectral  point 
will  be  approximately  25f.  Obviously  then,  as  N  is  reduced, 
the  spectral  points  become  more  spread  out  over  the  band, 
but  wider  in  bandwidth. 

Autocorrelators  can  also  be  implemented  using  analog 


or  digital  circuitry  or  combinations  of  the  two.  Figures 
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11  and  12  show  in  block  form  the  basic  functions  required 
for  each  implementation. 

3 . 3  Power  Spectrum  Estimate 

A  note  should  also  be  made  that  the  methods  of  spec¬ 
tral  analysis  described  here  do  not  provide  the  exact  power 
spectrum.  To  determine  the  average  power  of  a  stochastic 
signal,  an  infinite  time  would  be  required.  That  is,  T 
and  K  would  need  to  approach  infinity  in  Equations  (3)  and 
(9) .  Experimentally  obtaining  an  infinite  time  record 
is  not  a  reasonable  idea,  and  thus,  a  truncated,  or  esti¬ 
mate,  version  of  the  power  spectrum  is  accepted.  A  measure 
of  uncertainity  in  the  spectral  measurement  is  the  variance 
or  sensitivity,  and  is  shown  in  Equation  (1). 

Again,  the  idea  here  in  presenting  the  different 
forms  of  spectral  measurement  is  not  to  provide  a  vast 
theoretical  background.  A  rigorous  theoretical  development 
can  be  found  in  many  of  the  references  listed,  and  in 
particular,  from  a  most  basic  form  in  Lewis  (1977) . 

Instead,  the  intent  was  to  show  the  relative  complexities 
in  the  various  implementations  of  a  spectral  analysis. 
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CHAPTER  IV 

CHOICE  OF  SPECTRAL  ESTIMATION  TECHNIQUE 

4 . 1  Stability  Considerations 

The  basic  considerations  in  choosing  a  spectral 
estimation  technique  are  stability,  complexity,  and  cost. 
Stability  becomes  a  problem  in  analog  circuitry.  Amplifier 
and  multiplier  gains  must  be  kept  constant.  Filter  center 
frequencies  and  bandwidths  must  be  kept  constant.  Also 
integrator  time  constants  and  delay  line  delay  times 
must  remain  constant.  These  variables  must  not  only  re¬ 
main  constant  over  integration  times  on  the  order  of  hours, 
but  must  also  remain  equivalent  from  channel  to  channel  for 
perhaps  tens  of  channels.  The  major  stability  concern  for 
the  digital  circuitry  stems  from  keeping  the  sampling  inter¬ 
val  constant.  This  requires  only  a  single  oscillator  to 
remain  stable  over  time. 

4 . 2  Bandwidth  Considerations 

The  20  MHz  bandwidth  requirement  makes  the  digital 
circuitry  difficult.  High  speed  circuitry  must  be  used, 
thus  timing  becomes  critical.  Also  gate  delays,  crosstalk, 
and  even  transmission  line  reflections  become  important. 

The  digital  filter  method  has  an  additional  difficulty, 
in  that  to  attain  proper  channel  centers  and  bandwidths, 
a  significant  number  of  bits  are  required.  Therefore,  at 
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least  one  multibit  multiplier  is  required  per  channel, 
and  at  a  Nyquist  rate  of  40  MHz,  the  multiplication  must 
be  performed  in  less  than  25  nanoseconds.  The  digital 
autocorrelator  technique  alleviates  this  problem  since 
the  autocorrelation  function  can  be  performed  using  as 
few  as  one  bit.  Multiplications  now  became  only  an 
"exclusive  or"  function  and  can  be  performed  in  one  gate 
delay.  This  gross  quantization  error  is  paid  for  by  a 
loss  of  sensitivity.  The  correlation  time  must  then  be 
lengthened  to  achieve  a  sensitivity  equal  to  a  multibit 
autocorrelator  or  a  filter  bank. 

A  20  MHz  bandwidth  would  not  be  a  difficult  pro¬ 
cedure  for  analog  circuitry;  however,  the  stability  require¬ 
ments  would  cause  an  increase  in  the  circuit  complexity 
and  cost. 

4. 3  Cost  Cons ideration 

A  good  indication  of  relative  cost  is  found  by 
examining  the  "front  ends"  of  each  spectral  estimation 
technique.  The  final  accumulations  in  all  cases  can  most 
easily  be  done  using  a  digital  computer.  Therefore,  most 
of  the  expense  and  design  effort  would  be  expended  on  the 
high  speed  front  end  of  each  technique. 

A  typical  block  diagram  of  a  recursive  digital 
filter  is  summarized  in  Figure  13.  The  circuit  is  imple¬ 
menting  an  8  bit  version  of  the  difference  equation: 
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4 

y(n)  =  Z  y(n-i)a.+x(n-i)b. 

i=l  L  1 

where  the  input  data  stream  x(n-i)  is  quantized  only  to 
one  bit.  The  basic  cost  is  then  approximately  $240.00  per 
spectral  point  (Figure  13)  --  not  considering  the  control 
circuitry  required. 

A  one  by  one  bit  digital  autocorrelator  is  blocked 
out  in  Figure  14.  One  point  on  the  autocorrelation  function 
can  then  be  built  for  approximately  $30.00,  less  the  control 
circuitry  required. 

As  far  as  analog  circuitry  is  concerned,  crystal 
filters  sell  for  about  $250.00  each,  and  20  nanosecond 
delay  lines  are  between  $15.00  and  $40.00  each,  depending 
on  rise  times  and  insertion  losses. 

4 . 4  Digital  Autocorrelator 

Considering  the  system  stability  requirements  and  the 
complexity  and  cost,  a  digital  autocorrelator  appears  to 
be  the  most  attractive  spectral  estimation  technique.  This 
choice  restricts  the  spectral  point  parameters  to  being 
relatively  interdependent  ;  however,  uniformity  over  the  band 
is  then  easily  maintained.  Another  feature  of  the  system 
is  the  ease  of  changing  the  filters'  center  frequencies. 

A  simple  adjustment  of  the  sampling  rate  will  vary  At  and 
thus  effect  a  system  "fine  tuning."  This  is  attractive  for 
locating  the  peak  of  the  water  vapor  line. 
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Digital  autocorrelation  techniques  will  now  be  con¬ 
sidered  in  depth  and  a  particular  system  will  be  chosen. 


CHAPTER  V 


DIGITAL  AUTOCORRELATOR  TECHNIQUES 

5 . 1  Available  Equipment 

The  simplest  method  of  digital  signal  analysis 
would  be  a  complete  software  package.  The  Ionosphere 
Research  Laboratory  at  The  Pennsylvania  State  University 
has  available  a  Data  General  Nova  1200  minicomputer 
complete  with  an  IBM  nine  track  tape  interface.  In 
addition,  an  IBM  370-168  is  available  for  batch  processing 
at  Penn  State.  The  desired  system  would  be  to  digitize 
the  incoming  data  and  store  it  on  nine  track  tape  via  the 
Nova.  Data  analysis  could  then  be  completed  on  the  IBM  370. 
However,  the  20  MHz  IF  bandwidth  makes  a  hardware  accumulator/ 
buffer  necessary,  before  the  Nova  could  accept  the  data  stream. 
The  hardware  accumulator/buffer  is  to  be  a  high  speed  digital 
system  performing  the  autocorrelation  function  given  in 
Equation  (6) . 

5 . 2  Sampling  Rate 

The  sampling  rate  necessary  to  achieve  a  20  MHz  band 
would  be  40  MHz,  by  the  sampling  theorem.  However,  the 
sampling  theorem  assumes  an  infinite  number  of  quantiza¬ 
tion  levels.  Van  Vleck  (1943)  has  shown  that  the  auto¬ 
correlation  of  a  grossly  quantized  Gaussian  signal  can  be 
corrected  to  give  the  actual  autocorrelation  function. 


This  means  hardware  can  be  simplified.  However,  the  gross 
quantization  causes  a  decrease  in  sensitivity.  Hagan  and 
Farley  (1973)  and  Cooper  (1970)  show  that  an  improvement 
in  sensitivity  can  be  attained  by  over  sampling.  It  will 
be  shown  later  that  sampling  at  twice  the  Nyquist  rate , 

80  MHz,  gives  nearly  maximum  improvement  in  sensitivity. 

Also,  if  a  frequency  "fine  tuning"  is  desired,  the  sampling 
rate  should  be  further  increased  to  allow  At  to  vary  with¬ 
out  decrementing  system  sensitivity.  This  concept  will  be 
discussed  further,  but  for  now  it  is  sufficient  to  assume  that 
the  system  should  be  designed  for  a  sampling  rate  of  100  MHz. 

5 . 3  Number  of  Quantization  Levels 

The  next  step  is  to  consider  the  number  of  quanti¬ 
zation  levels  desired.  A  large  number  of  bits  would  de¬ 
crease  the  sampling  frequency  down  to  around  40  MHz,  how¬ 
ever,  even  then  a  multiplication  operation  must  be  accom¬ 
plished  in  less  than  25  nsec.  P.  Blankenship  and  A.  Huntoon 
have  done  some  digital  processing  work  (Rabiner  and  Gold. 

1975)  and  have  accomplished  a  9  by  9  bit  multiply  in  22 
nsec  using  45  custom  2  bit  ECL  adders.  On  the  other  hand, 
a  number  of  1  by  1  bit  autocorrelators  are  presently  in 
operation  (Weinreb,  1963;  Fell  and  Little,  1975;  Abies, 

Cooper,  Hunt,  Moorey,  and  Brooks,  1974).  The  1  bit  quanti¬ 
zation  requires  a  multiply  every  10  nsec,  however,  a  lxl 
bit  multiplication  is  simply  an  "exclusive  or"  function 
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and  can  be  accomplished  by  the  ECL  III  logic  series  in 
one  nsec.  Therefore,  a  tradeoff  exists  between  hard¬ 
ware  complexity  and  sensitivity.  An  interesting  compromise 
is  a  1  x  2  bit  system  where  the  one  bit  quantization  is  also 
the  delayed  signal.  The  1x2  bit  system  still  maintains 
a  sign  only,  or  one  bit,  multiplication.  Also,  the 
delayed  signal  being  one  bit  means  only  one  channel  of  shift 
registers  is  required.  The  1x2  bit  system  maintains  much 
of  the  simplicity  offered  by  the  lxl  bit  system,  however, 
the  sensitivity  can  approach  a  maximum  of  847„  of  a  multibit 
correlator  (Hagan  and  Farley,  1973),  as  opposed  to  74% 
(Hagan,  1972)  in  the  lxl  case.  The  1x2  bit  digital 
autocorrelator  was  then  accepted  as  the  overall  best  system 
for  analyzing  the  water  vapor  line. 


aA 


—  m 
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CHAPTER  VI 

ONE  BY  TWO  BIT  CORRELATOR 

6 . I  Normalized  Autocorrelation  Function 

The  1x2  bit  autocorrelation  estimate,  as  Van  Vleck 
showed,  contains  a  gross  quantization  error  and  therefore 
must  be  corrected  to  provide  an  equivalent  multibit  auto¬ 
correlation  estimate.  This  correction  factor  can  be 
arrived  at  by  considering  the  statistical  mean  of  the 
normalized  autocorrelation  function  (uncorrected) ,  r(iAt) . 
Note  that  the  normalized  autocorrelation  function  is 
referenced  to  the  zero  lag  autocorrelation  point. 

r(il,)  -  5351 

From  Equation  (6) : 

r(iAt)  =  S-feAt)  -x  OjAt+iAx)  (13) 

x2 (kAt) 

Here  x(kAt)  represents  the  two  bit  real  time  sig¬ 
nal  and  x(kAt+iAT) ,  represents  the  one  bit  delayed  signal. 
The  mean  of  the  signal  products  can  then  be  represented 
by  the  sum  of  the  joint  probabilities  of  each  quantization 
product  occurring.  The  quantization  products  for  1x2 
bit  quantization  scheme  are  defined  in  Table  1.  Note  that 
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Table  1.  One  by  Two  Bit  Quantization  Products 

States : 

One  bit  [x(kAt  +  iAr) ] 

1  x  >  0 

1  x  <  0 


Two  bit  [x(kAt) J 

2  V  <  x  <  » 

o 

1  0  <  x  <  V 

o 

1  -V  <  x  <  0 

o 

2  -oo  <  x  <  -V 

o 


Products : 

Two  bit 


One 

bit  _ 

1  -1  -n  n"1  1 
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the  quantization  LeveL,  v  ,  is  normalized  to  the  rms  value 
of  the  incoming  signal.  Also  the  higher  products  are 
weighted  by  n.  The  proper  choice  of  n  will  tend  to  im¬ 
prove  sensitivity. 

The  probability  of  a  particular  product  occurring 
is  then  determined  by  the  joint  probability  of  the  corres¬ 
ponding  delayed  and  real  time  states  occurring.  Each 
probability  function,  Pj  ,  is  denoted  by  the  subscripts  d 
and  r,  indicating  the  state  of  the  delayed  and  real  time 
signals,  respectively.  Equation  (13)  can  then  be  repre¬ 
sented  as  follows: 

r(iAt)-t  (l)  (pl2+pi:>  +  (-D  (P12+Pl2)  +  (n"1)  (P11+Pll)  + 

C-n”1)  <PL1+Pll>  1  '  t  CD  CPL2+P-5)-H<n'1) 

<pu+pn>> 

Assuming  the  input  signal  is  a  Gaussian  random  process  with 
zero  mean,  symmetry  will  allow  the  following  reduction. 


(Pi2-pi2)+Tl  (pirpii) 

r(iAx)  =  — it — - — — — 


(14) 


P12+n  P11 


The  joint  probability  function  is  then  defined  by  the 
Bivariate  Normal  Probability  Function  (National  Bureau  of 
Standards,  1964)  where  p  is  the  normalized  correlation. 
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v  7  (  n  ^  (y\  VJ.1 

Par  -  QU)Q(r)^o  | 

(15) 

Where : 

Q(x) 

r  Z(c)dt 

*  X 

(16) 

Z(x) 

r  In 

(17) 

Z(m)(x)  - 

.m 

£=zw 

(18) 

Some  supporting 

probability  functions  are  also 

defined : 

P(x) 

rx 

Z(t)dt 

J  -  oo 

(19) 

A(x) 

f  Z(t)dt 

J  -X 

(20) 

P(x)+Q(x) 

-  1 

(21) 

P(-x) 

-  Q(X) 

(32) 

A(x) 

-  2P(x) -1 

(23) 

Combining  Equations  (31)  and  (33)  then  provides 

: 

Q(x) 

-  j  -  j  A(x) 

(24) 

Where  A(x)  has  been  tabulated  for  various  values  of  x  in 
a  variety  of  handbooks  (Dwight,  1947). 


The  probability  functions  P^r  can  now  be  determined 
to  any  accuracy  by  setting  the  limit  on  N  ir.  Equation  (15)  . 
For  the  probability  P,.-,,  Equation  (15)  is  used  in  conjunc¬ 
tion  with  Equations  (16)  through  (24).  Here  N  is  set  to  4. 


1  Z(V  ) 

P12  “  i(l-A(Vo)-K_^_p 


6.2? 


[Z(Vo)](Vo4-6Vo3+3)  5 


40. 


(25) 


Obviously  as  p  decreases,  that  is,  the  signals  become  less 
correlated,  the  higher  order  terns  become  less  significant. 

The  probability  P-j-.,  has  limits  -*>  to  0  for  the  d 
parameter.  Therefore,  Equation  (22)  can  be  substituted 
into  Equation  (15)  giving: 


> 


Pdr  =  ?(d>Q(r)+ 


,(N) 


-d): 


N=0 


(N) 


N-l)  ! 


iEl  ,*-1 


(26) 


Then  substituting  as  before  gives: 

l  Z(V  )  Z(V)(V  :-l)  3 

Pl2  =  ^l-A(Vo)  1-— £-o  +  - -  p3 


/2i 

Z(VoHVo4-6V03+3)  5 

40/2? 


0.  2i 


(27) 


P^^  can  be  found  by  subtracting  from  the 

probability  defined  by  both  limits  dir  going  from  0  to 
*.  These  relationships  are  then: 


00  £ 


+  27P  +  L  -3 


nrp' 


+  2^37>v 5 


therefore : 
P 


11 


1  l  Z(V  ) 

5a<Vo>+<^  -  -=^>0  + 


v  ZTT 


:(v  ) 


[yr  +  - — '  2 


r  Z7T 


(vo  -i)] 


P3 

5~ 


3  2<V 


+  ^  *  -^<V04-6V02+3)1&5 


✓27 


is 


11  ls  approached  similarly: 


P11  ^O^U 


poo 


1 

4 


1  13  g 

"  r27p  -  2£S7P 


and  therefore: 


JA<Vc>-(2T 
3 


Z(V„) 


—  4-  f  -2 

6  +  27 


r  /  7T 

Z(VJ 


21)  +  — 


✓27 


✓27 


°^Vo4-6Vo2+3>]fe 


) 


The  relationship  between  the  normalized  uncorrected 
autocorrelation  function,  r(iAx),  and  the  corrected  version, 
p,  is  found  by  substituting  Equations  (25),  (27),  (28),  and 
(29)  into  Equation  (14) .  The  rather  forbidding  result 
can  then  be  solved  for  p  by  a  computer,  using  iterative 
methods.  However,  a  simple  approximation  would  be  to 
assume  highly  uncorrelated  signals,  as  would  be  the  case 


Substituting  Equations  (30)  through  (33)  into  Equation 
(14)  now  produces: 
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It  is  desirable  to  have  n  a  power  of  two  so  that  weightings 
can  be  easily  achieved  with  a  digital  counter.  With  Vq  and 
n  optimized,  integration  times  relative  to  a  multibit  corre¬ 
lator  were  found  to  be  1.78  for  sampling  at  the  Nyquist 
rate,  and  1.43  for  oversampling  at  twice  the  Nyquist  rate. 

A  sampling  rate  higher  than  twice  the  Nyquist  rate  does 
afford  some  further  increase  in  sensitivity  (Burns  and 
Yao,  1969;  Hagan,  1972).  However,  the  slight  sensitivity 
increase  does  not  justify  a  further  increase  in  bandwidth. 

It  is  adequate  to  realize  that  sampling  faster  than  two 
times  the  sampling  rate  does  not  degrade  sensitivity. 

The  relative  integration  times  determined  by  Hagen 
and  Farley  can  be  converted  to  relative  sensitivities 
using  Equation  (1) .  Since  sensitivity  is  inversely  pro¬ 
portional  to  the  square  root  of  integration  time,  sampling 
at  the  Nyquist  rate  will  provide  a  relative  sensitivity 
of  0.75,  and  twice  the  Nyquist  rate  will  increase  the 
sensitivity  to  .84  of  a  multibit  correlator. 

To  construct  the  most  sensitive  2x1  bit  correlator 
capable  of  analyzing  a  20  MHz  bandwidth,  a  minimum  sampling 
rate  of  two  times  the  Nyquist  rate  or  80  MHz  is  to  be  used. 
Then,  as  previously  discussed,  to  allow  some  flexibility  in 
the  location  of  the  resultant  "spectral  windows,"  the 
sampling  rate  will  be  made  variable  between  80  and  100  MHz. 
Oversampling  can  be  accomplished  by  setting  4t  equal  to 


45 


2At  in  Equation  (6) .  This  would  be  equivalent  to  delaying 
the  signal,  x(kAt  +  iAr) ,  through  two  successive  At  clocked 
flip  flops  before  i  is  considered  incremented.  Weighting 
the  higher  level  products  by  four  is  most  easily  accomplished 
by  considering  the  higher  level  products  one,  and  dividing 
the  lower  level  products  by  four.  An  up-down  counter  will 
then  accumulate  the  products.  A  simple  block  diagram  of 
this  system  is  represented  in  Figure  12. 
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CHAPTER  VII 
SYSTEM  HARDWARE 

7 . 1.  Logic  Families 

The  fastest  Schottky  TTL  up-down  counter  listed  in 
the  Texas  Instruments  TTL  Data  Book  will  clock  at  a  maximum 
rate  of  40  MHz.  This  makes  the  TTL  logic  series  inadequate 
for  the  high  speed  section  of  the  correlator.  The 
Motorola  ECL  Data  Book  lists  several  ECL  logic  series  which 
feature  much  faster  gate  delay  times  then  the  TTL  families. 
The  ECL  III  series  lists  a  1  nsec  typical  gate  delay  time 
as  opposed  to  the  3  nsec  delay  time  offered  by  the  Schottky 
clamped  TTL.  Motorola  also  offers  an  ECL-10,000  series 
with  gate  delay  times  of  2  nsec.  The  10,000  series  pro¬ 
vides  gate  rise  and  fall  times  intentionally  slowed  to 
2  nsec  to  help  inhibit  coupling  between  adjacent  signal 
lines.  The  MC  10136  up-down  counter  lists  a  150  MHz 
maximum  clocking  rate  in  the  up  or  down  mode.  If  the  up 
and  down  modes  are  used  alternately,  a  maximum  clocking 
rate  of  133  MHz  is  inferred.  Also,  cascading  the  counters 
slows  the  maximum  clocking  rate  to  95  MHz.  An  additional 
feature  of  the  ECL  10,000  series  is  the  variety  of  MSI 
and  LSI  functions.  The  readily  available  array  of  complex 
functions  is  comparable  to  the  Schottky  TTL  series.  Thus, 
the  ECL  10,000  series  appears  adequate  to  provide  the 
necessary  high  speed  logic.  The  higher  cost  associated 
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with  ECL  ($10.50  for  a  four  bit  binary  up-down  counter) 
will  then  restrict  ECL  to  only  the  high  speed  section. 
Subsequent  lower  speed  sections  will  then  be  built  from 
the  less  expensive  TTL  families. 

The  ECL  10,000  series  consists  of  open  emitter 
follower  output  stages.  Therefore,  with  the  proper  pull 
down  resistor,  each  logic  gate  can  become  a  line  driver  or 
each  gate  can  drive  low  impedance  transmission  lines,  down 
to  an  impedance  of  50  ohms.  At  switching  times  of  100 
MHz,  transmission  line  effects  become  important  and,  thus, 
should  be  considered  in  circuit  board  layouts. 

First,  however,  the  power  supply  voltage  should  be 
considered.  According  to  Motorola's  ECL  data  book,  a 
supply  voltage  of  -5.2  volts  will  provide  best  noise 
immunity.  The  resulting  logic  levels  are  then  -0.9  volts 
for  a  logic  high  or  "1"  and  -1.75  volts  for  a  logic  low 
or  "0."  The  resulting  logic  level  transition  voltage  is 
-1.29  volts. 

7 . 2  Crosstalk 

Crosstalk  should  be  considered  at  these  switching 
speeds  and  is  easily  handled  by  using  a  method  described 
in  an  article  by  Robert  Saeger  of  Signetics  Corporation 
(1974) .  Jaeger  describes  the  relationship  between  the 
sending  and  receiving  line  in  the  form  of  pulse  widths  and 
amplitudes  induced  onto  the  receiving  line.  For  the 
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purposes  of  the  system  here,  pulse  widths  are  fairly  unim¬ 
portant,  provided  the  induced  pulse  amplitudes  are  small. 
Crosstalk  is  broken  down  into  a  forward  going  pulse  rela¬ 
tive  to  the  sending  line  and  backward  pulse.  The  forward 
and  backward  crosstalk  amplitudes  were  determined  by  Jaeger 
to  be : 


Vf 


KPV 
f  s 

tr 


V.  = 


K.Vo 
b  s 


(37) 


(38) 


Where : 

Vf  =  Amplitude  of  forward  pulse 
Vg  =  Amplitude  of  sending  pulse 
=  Amplitude  of  backward  pulse 
i  =  Length  of  parallel  sending  and  receiving 
lines  in  feet 

tf  *  Sending  pulse  rise  time  in  nsec 
Kg  =  Forward  crosstalk  coefficient 
K^  =  Backward  crosstalk  coefficient 


The  coefficients  Kg  and  K^  are  dependent  on  the  spacing 
between  sending  and  receiving  lines  and  are  provided  in 
chart  form  by  Jaeger  and  repeated  here  in  Figure  15a  and  b. 


Now,  by  first  considering  the  backward  crosstalk, 
a  minimum  line  spacing  of  100  mils  provides  about  the  most 
effective  crosstalk  coefficient  while  still  providing  a 
reasonable  line  spacing.  From  Equation  (37)  ,  a  of 
0.05  and  a  Vg  or  logic  level  difference  of  0.85  volts 
will  then  give  a  of  0.043  volts  which  is  insufficient 
to  cause  either  logic  level  to  approach  the  logic  cross¬ 
over  voltage. 

Then  considering  forward  crosstalk,  100  mil  line 
spacing  gives  Kf  as  0.055.  For  ECL  10,000  logic,  t  is 
2  nsec  and  assuming  a  crosscoupling  length  of  1  ft,  Vf 
becomes  0.023  volts.  Forward  crosstalk  is  therefore  less 
of  a  problem  then  backward  crosstalk. 


7 . 3  Transmission  Lines 

Another  problem  associated  with  fast  switching 
lines  is  ringing  or  reflections  on  a  signal  line  due  to 
impedance  mismatching.  From  transmission  line  theory, 
the  reflection  coefficient  is  given  as-. 

KR  =  ZV+Z° 


Where : 

=  Reflection  coefficient 
ZQ  =  Transmission  line  impedance 
=  Termination  or  load  impedance 


The  easiest  mathod  to  inhibit  ringing  then  is  to  make 
Z^  =  ZQ.  Thus,  the  value  of  the  pull  down  resistors  for 
each  logic  output  is  dependent  on  the  signal  line  impedance. 
The  value  of  ZQ  can  be  estimated  from  a  microstrip  repre¬ 
sentation  of  the  circuit  board.  Using  a  glass-epoxy  circuit 
board  of  width  0.0625"  and  permittivity  of  4.6,  the  impedance 
can  be  determined  from  Figure  16.  The  impedance  is  depedent 
on  signal  line  width,  assuming  the  line  thickness  is  small 
compared  to  width.  A  typical  line  thickness  will  be  approx¬ 
imately  1.5  mils,  whereas  the  smallest  easily  etchable 
widths  will  be  10  to  20  mils. 

Since  the  ECL  10,000  series  was  intended  to  drive  low 
impedance  loads,  a  wide  line  width  is  desired.  The  widest 
pratical  line  width  for  connections  to  DIP  circuit  packages 
is  50  mils.  If  a  ground  plane  could  be  managed,  a  ZQ  of 
75  ohms  would  be  maintained.  However,  two  sided  circuit 
boards  were  required  in  the  circuit  layouts,  therefore  only 
partial  ground  grid  patterns  were  constructed  on  both  sides 
of  the  board.  The  resulting  effect  would  be  to  increase 
the  average  line  impedance  and  lower  the  uniformity  along 
the  line.  This  lack  of  carefully  controlled  line  impedance 
may  seem  crude,  however,  Jaeger  suggests  that  the  EC1.  10,000 
series  will  work  even  if  wire-wrapped  connections  are  used. 

The  actual  terminations,  or  pull  down  resistors  were 
found  experimentally.  A  Tektronics  475  oscilloscope  was 
the  best  oscilloscope  available,  however,  the  bandwidth 
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was  only  200  MHz.  This  made  the  detection  of  ringing  on 
a  100  MHz  square  wave  impossible.  Thus,  the  terminations 
were  made  220  Q,  which  gave  the  highest  amplitude  signal 
at  the  receiving  ends  of  a  majority  of  transmission  lines. 

One  other  circuitry  problem  is  the  connections 
between  circuit  boards.  Because  of  the  lack  of  good  ground 
planes,  chip  to  chip  connections  between  boards  were  made 
for  high  speed  signals  using  twisted  pair  lines.  Number 
24  wire  with  approximately  13  turns  per  foot  was  used. 

The  following  chapters  will  describe  the  actual 
circuitry  used. 
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CHAPTER  VIII 
CIRCUITRY 


8 . 1  Quantizer 

The  quantizer  was  constructed  on  the  timing  genera¬ 
tion  board  (TG)  utilizing  two  Motorola  MC  1650  comparators. 
These  comparators  have  a  compare  time  of  3.5  nsec  and 
include  D  flip-flop  outputs.  Although  each  ECL  10,000 
series  output  is  capable  of  driving  up  to  90  gates,  rise 
time  does  increase  with  fanout.  Thus,  fanout  was  kept  to 
a  maximum  of  4,  and  the  sytem  was  designed  to  accommodate 
16  autocorrelation  points.  This  circuitry  is  shown  in 
Figure  17. 

Note  here  that  V  and  -V  were  set  to  +.5  v.  and 
o  o 

-.5  v.,  respectively.  This  indicates  an  analog  signal  of 

475  mv.  rms  is  required  to  provide  the  system  sensitivity 

previously  discussed.  The  accuracy  of  V  will  affect  sys- 

o 

tem  sensitivity  according  to  Equation  (36) .  From  the  work 
of  Hagan  and  Farley  (1973) ,  Equation  (36)  will  result  in  an 
integration  time  relative  to  a  multibit  correlator  and  is 
plotted  in  Figure  18  verses  Vq  superimposed  in  the  relative 
sensitivity  as  derived  from  Equation  (1) .  Obviously  the 
value  of  Vq  is  not  critical.  A  10%  change  in  the  optimum 
VQ  will  result  in  only  a  0.13%  degradation  in  sensitivity. 

8 . 2  High  Speed  Accumulator 

The  high  speed  section  (HS)  runs  synchronously 
from  a  100  MHz  master  clock.  Figure  19  shows  the  delay 
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line  and  multiplier.  The  delay  line,  an  MC  10141  shift 
register,  will  allow  two  autocorrelation  points  if  over- 
sampling  is  used.  Thus,  in  Equation  (6),  At  =  2At  =  20  nsec, 
and  therefore  two  points  were  laid  out  on  each  board.  Re¬ 
lating  to  the  labeling  of  circuit  diagrams,  if  two  integrat¬ 
ed  circuit  number  appear,  the  upper  number  refers  to  point 
i  and  the  lower  number  to  i  +  1.  The  i=0  and  i=l  boards 
are  special  in  that  only  one  delay  is  required.  This 
special  wiring  is  indicated  by  the  dotted  lines  and  break 
marks.  If  operation  at  only  the  Nyquist  rate  is  preferred, 
the  connections  to  pins  15  and  3  on  MC  10141  can  be  broken 
and  jumped  to  14  and  15,  respectively.  This  will  allow 
operation  with  At  •  At. 

The  delayed  and  real  time  signal  multiplication  is 
only  a  sign  multiplication  and  is  accomplished  using  10107 
"exclusive  or"  gates.  The  following  D  flip-flops  serve 
only  to  retain  synchronous  operation  before  the  resultant 
sign  and  magnitude  bits  are  combined  to  produce  high  and 
low  level  counts,  both  up  and  down,  as  per  Table  1.  Note 
that  the  higher  level  products  being  weighted  by  n  (recall 
n=4)  is  accomplished  by  dividing  the  lower  level  products 
by  n.  Figure  20  shows  the  divide  by  4  circuitry  built  from 
J/K  flip-flops.  The  divide  by  4  counter  was  constructed  as 
a  Mealy  machine  and  is  described  in  Figure  21  by  means  of 
the  flow  diagram  and  coded  state  table.  The  outputs  were 
then  or'ed  with  the  high  level  counts  and  the  result 
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K,  = 

xy2  .  . 

Kr  X  +  Q2 
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k2* 

x 

Kf* 

Z  = 

YYX  .  . 

Z  =  Q  +  Q  +  X 

figure  21.  Divide  by  4  counter  Logic 


manipulated  to  fit  an  MC  10118  bilevel  gate.  The  counts 
are  then  accumulated  in  an  MC  10136  binary  up-down  counter 
shown  in  Figure  22. 

Because  of  the  high  expense  of  ECL  counters,  a 
transfer  was  made  to  TTL  at  this  point.  The  MC  10136  is 
operated  as  a  divide  by  16  counter  carrying  only  one  bit 
to  the  lower  speed  accumulators.  At  the  start  of  each  16 
count  cycle,  the  counter  is  preloaded  to  8  and  then  tested 
for  a  positive  or  negative  pass  through  zero.  A  positive 
or  negative  count  greater  than  eight  is  then  carried  to 
the  TTL  counters,  facilitating  a  rounding  operation.  The 
reset  pulse  which  then  reloads  the  counter  to  eight  is 
entered  via  the  output  flip-flops  shown  in  Figure  20. 

The  zero  detection,  denoted  x,  had  to  be  derived  from  the 
counter  outputs  since  the  carry  bit  used  in  cascading  has 
too  great  a  propagation  delay  time.  The  present  count's 
sign,  denoted  s,  was  easily  derived  from  the  up-down  inputs 
with  a  single  J/K  flip-flop.  The  logic  for  determining  a 
TTL  count  and  its  corresponding  sign  is  shown  in  Figure  23 
in  the  form  of  a  flow  diagram  and  coded  state  table. 

The  effect  of  using  what  amounts  to  a  second  quanti¬ 
zation  scheme  here  can  intuitively  contribute  to  a  further 
decline  in  sensitivity.  In  actuality,  however,  the  decline 
in  sensitivity  can  be  made  negligible  if  the  incoming  signal 
is  somewhat  correlated  and  the  number  of  samples  is  very 
large.  Appendix  A  provides  further  analysis  of  this  type  of 
rounding  error. 
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Figure  23. 

High  speed  counter  logic 

8 . 3  Signal  Propagation 

The  use  of  two  sided  boards  as  opposed  to  three 
sided  (center  ground  plane)  in  the  high  speed  section 
caused  some  propagation  problems,  especially  of  the  high 
frequency  clock  signal.  Board  to  board  connections  were 
done  simply  by  carrying  a  ground  from  point  to  point  along 
with  the  signal  in  a  twisted  pair  line.  On  board  propa¬ 
gation  of  50  MHz  or  lower  signals  (non-clock)  proved 
not  to  be  a  problem  when  terminated  at  a  single  destination 
with  220  0.  The  100  MHz  clock  signal  however,  proved  more 
difficult  in  some  instances  which  were  rectified  as  follows 

1.  Occasionally  reflections  at  the  source  only 
will  be  a  problem,  thus  a  source  termination  of  220  0  is 
effective . 

2.  When  reflections  at  both  source  and  destination 
are  a  problem  or  crosstalk  signals  appear,  a  destination 
terminated  with  220  0  and  the  source  terminated  with  a 
small  capacitance  (10-20  pf)  is  effective. 

3.  For  badly  distorted  signals,  a  twisted  pair  line 
directly  connecting  source  to  destination  and  terminating 
the  destination  in  200  0  transmits  the  signal. 

4.  If  long  board  to  board  connections  cause  ex¬ 
cessive  attenuation,  a  differential  twisted  pair  can  be  run 
using  ECL  chips  10115  and  10101. 
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When  many  high  speed  boards  are  used,  approximately 
equal  length  lines  are  required  to  transmit  the  various 
signals  to  each  board.  Proper  operation  of  each  board  can 
be  determined  by  placing  a  short  pulse  (10  nsec)  on  the 
analog  input  then  measuring  the  time  delay  between  the 
real  time  and  delayed  sign  bits.  This  measurement  can  be 
made  at  HS  13  pin  7  and  HS  13  pin  9  or  their  equivalents. 
The  delayed  pulse  should  follow  the  real  time  pulse  by  the 
proper  number  of  delays.  Adjustment  of  the  delayed  time 
pulse  can  be  made  by  "phasing”  (delaying)  the  clock  signal 
at  that  particular  board  using  one  or  more  gates.  These 
gates  are  provided  on  HS  14a.  The  real  time  signal  should 
not  be  delayed  as  this  signal  will  be  in  phase  with  the 
magnitude,  reset,  and  enabling  pulses. 

8 . 4  Medium  and  Low  Speed  Accumulators 

Four  medium  speed  accumulators  with  twelve  bit  tri¬ 
state  latches  are  included  on  one  printed  circuit  board. 
Figure  24  shows  the  circuitry  with  the  package  and  pin 
numbers  shown  in  quadruplicate.  The  top  number  refers  to 
channel  i  and  then  in  descending  order  to  channel  i  +  3 
as  referenced  to  Equation  (6) .  These  74LS191  up-down 
counters  accumulate  the  counts  produced  by  the  high  speed 
section  and  are  synchronously  clocked  by  the  100  MHz  master 
divided  by  18.  Although  only  16  counts  were  made  in  the 
high  speed  section,  two  counts  were  required  to  allow  for 
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the  reset  pulse,  thus  forcing  the  medium  speed  counters  to 
be  clocked  at  1/13  speed.  This  will  also  decrease  sensi¬ 
tivity  by  /1/9  since  1/9  of  the  samples  are  not  used. 

The  12  bit  word  accumulated  is  in  two's  compliment  form, 
the  most  significant  bit  being  a  sign  bit.  Therefore, 
the  counters  are  operated  as  divide  by  2.  loading  the 
contents  into  three  74173,  four  bit  latches  every  3.6864 
nsec.  The  latch  outputs  are  then  all  bussed  together  to 
form  the  12  bit  data  link  (DL4>  to  DL11)  . 

The  low  speed  accumulator  is  built  on  board  one  of 
two  wire  wrap  boards.  This  accumulator  uses  eight  Intel 
2111A-4  RAM's  producing  256  32-bit  words.  A  256  deep 
RAM  allows  for  an  eventual  system  expansion  to  256  points. 
With  the  RAM  cycle  time  of  only  450  nsec,  this  not  only 
allows  all  RAM  words  to  be  accessed  but  also  allows  approx¬ 
imately  13.9  ps  to  manipulate  each  word  before  the  medium 
speed  latches  are  reloaded.  The  32-bit  word  width, 
accessed  via  a  74LS257A,  2  to  1  multiplexor  allows  most 
efficient  access  from  16  or  8  bit  computer  systems.  Also, 
using  a  two's  complement  word  configuration,  the  32-bits 
will  accumulate  a  maximum  of  1048576,  12  bit  two’s  complement 
words.  In  accumulation  time,  this  amounts  to  about  64 
minutes.  Thus,  interfacing  to  most  computer  systems  should 
not  be  a  difficult  procedure. 
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The  RAM  word  is  added  to  the  incoming  12-bit  word  as 
shown  in  Figure  25.  The  12-bit  word  can  be  added  or  sub¬ 
tracted  from  the  running  sum  as  desired  via  the  invert  con¬ 
trol.  A  low  on  this  line  adds  while  a  high  subtracts.  The 
invert  control  is  also  connected  to  the  least  significant 
carry  input  on  the  32-bit  full  adder.  This  facilitates 
adding  one  after  a  complement  to  provide  a  two's  complement 
negation.  Thus,  the  RAM  word  is  read,  added  to  the  12-bit 
data  bus,  and  the  result  stored  in  the  32-bit  data  bus. 

The  32-bit  2  to  1  multiplexor  then  allows  the  sum  to  be 
writcen  back  into  RAM,  or  allows  the  RAM  to  be  cleared  by 
reloading  with  the  present  12-bit  word. 

The  computer  interface  logic,  also  shown  in  Figure 
25,  is  built  on  board  two  of  the  wire-wrap  set.  Six  74365 
tri-state  drivers  form  a  bus  receiver  which  isolates  the 
computer  action  from  the  continuous  correlator  operation. 
Before  each  word  in  the  low  speed  accumulator  (RAM) , 
discussed  above,  is  reloaded,  the  accumulator  read  phase 
writes  the  existing  sum  into  a  second  32  x  256  buffer  RAM. 
The  select  pulse  controls  this  operation  and  allows  the 
accumulator  RAM  to  be  loaded  into  the  buffer  RAM  at  least 
once  every  64  minutes,  as  previously  discussed.  Removing 
the  select  pulse  then  turns  control  of  the  buffer  RAM  over 
to  the  computer  through  two  sets  of  multiplexors.  The 
present  system  is  designed  to  be  accessed  from  a  16  bit 
Nova  1200  built  by  Data  General  Corporation.  Using  the 
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Figure  25.  Low  speed  accumulator  and  interface 


Nova's  least  significant  address  bit  as  a  selector,  the 
upper  and  lower  16  bits  are  accessed  alternately  through  a 
2  to  1  multiplexor.  The  Nova  can  then  further  accumulate 
for  the  experiment  duration. 

The  Nova  is  also  interfaced  to  an  IBM  compatible  9 
track  tape  drive.  Thus,  at  the  conclusion  of  the  experiment, 
data  can  be  transferred  to  an  IBM  370  computer  in  a  totally 
or  partially  accumulated  form.  The  final  accumulation,  if 
necessary,  along  with  the  Fourier  Transform  can  then  be 
accomplished  on  the  370. 

8 . 5  Timing 

The  timing  section  generates  the  wave  forms  required 
to  execute  the  loading  and  resetting  used  in  the  medium  and 
high  speed  sections.  Also,  the  waveforms  used  to  control 
the  RAM's  in  the  low  speed  and  interface  sections  are  gen¬ 
erated.  This  circuitry  is  built  on  the  timing  generation 
(TG)  board,  fanout  boards  one  and  two  (TFI,  TFII) ,  and 
board  two  of  the  wire  wrap  pair. 

The  reset  pulse  used  to  reload  the  high  speed  counter, 
shown  in  Figure  22,  actually  takes  two  10  nsec  clock  periods. 
Thus,  16  pulses  are  counted  and  two  used  to  reset.  A  divide 
by  18  counter  is  then  required  to  generate  the  reset  pulse 
as  shown  in  Figures  26a  and  26b.  Here,  the  10136  counter 
is  forced  to  stop  counting  for  two  pulses  via  the  associated 
flip-flops  and  gates.  The  flow  diagram  and  coded  state 
table  describing  this  circuitry  is  shown  in  Figure  27. 


FLOW  DIAGRAM* 
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Figure  27. 
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The  two  most  significant  counter  bits  are  then 
level  shifted  to  TTL.  Q2TTL  then  clocks  the  medium 
speed  counters .  Although  Q2TTL  has  two  pulses  for  every 
18  high  speed  counts  (see  timing  diagrams  in  Appendix  B) , 
one  pulse  occurs  during  the  high  speed  count  period  when 
the  medium  speed  counters  are  not  enabled.  Clock  Bl  is 
generated  from  Q2TTL  and  Q3TTL  (see  Appendix  B) . 

Clearing  the  medium  speed  counters  and  setting  the 
12  bit  latches  must  occur  before  the  medium  speed  counters 
overrun  the  sign  bit.  This  allows  a  maximum  of  2047 
counts  before  the  low  going  load  pulse  loads  the  present 
count  into  the  74173  latches.  The  load  occurs  on  the 
load  pulse  rising  edge  which  also  triggers  the  low  going 
clear  pulse.  The  clear  pulse  low  level  then  clears  the 
74LS191  counters.  Figure  28  shows  the  divide  by  2047 
counter  used  to  produce  the  clear  and  load  pulses. 

Clock  Bl  and  the  stretched  version,  Clock  B2,  are  derived 
from  Q3TTL  and  Q2TTL  to  clock  the  divide  by  2047  counter. 
Clock  Bl  had  to  be  stretched  in  order  to  accommodate  the 
the  short  ripple  output  of  TG22.  Note  at  this  stage  an 
external  reset  line  is  provided.  A  low  on  this  line  will 
allow  the  high  speed  section  to  fill  with  new  information 
when  the  analog  source  is  switched  to  a  reference.  The 
delay  involved  to  reload  the  high  speed  section  will  be 
less  than  300  nsec  for  16  points. 
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In  addition  to  the  clear  and  load  pulses,  the  divide 
by  2047  counter  generates  the  waveforms  used  to  control 
the  low  speed  accumulator  and  buffer  RAM.  Also,  the  address 
bus  is  generated  which  accesses  the  two  RAMs  and  is  decoded 
to  provide  access  onto  the  12  bit  data  link.  This  logic 
is  shown  in  Figure  28  with  the  accompanying  waveforms  given 
in  Appendix  B.  Figure  29  shows  the  address  decode  which, 
in  turn,  will  enable  each  medium  speed  latch  to  access  the 
12  bit  data  link. 

Figure  30  shows  the  select  pulse  generation,  used 
in  the  low  speed  accumulator  and  interface  sections. 

Here  load  pulses  are  counted  and  can  be  accumulated  up  to 
1048575  pulses  before  a  select  pulse  is  generated.  The 
select  pulse  will  last  through  one  cycle  of  the  medium 
speed  counter  which  amounts  to  368.38  nsec .  Presently 
the  system  is  wired  to  provide  a  select  every  262143  counts 
or  about  every  96.5  seconds.  This  allows  over  a  minute  for 
the  Nova  1200  to  access  16  locations  in  the  buffer  RAM. 

The  end  of  the  select  pulse  will  then  serve  to  interrupt 
the  Nova,  instructing  the  computer  to  begin  addressing 
the  buffer  RAM.  In  addition,  a  manual  high  on  the  start 
line  here  will  clear  this  counter  and  allow  the  initial 
11/2  minute  count  frame  to  start. 

If  other  count  periods  are  desired,  chips  67,  69, 
and  76  can  be  rewired.  Using  output  Q16  on  chip  76  will 
provide  a  cycle  time  of  24  seconds,  Q17:  48  seconds, 


Figure  29.  RAM  address  decoding 
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Figure  30.  Cycle  timing 
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Q18 :  96.5  seconds,  Q19 :  193  seconds,  and  using  no 

outputs,  that  is,  allowing  the  counters  to  completely 
cycle  through  zero,  will  provide  386  seconds.  Appendix  B 
gives  the  timing  diagrams  describing  the  entire  system. 
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CHAPTER  IX 
SYSTEM  OPERATION 

9 . 1  Initialization 

The  initial  start  up  of  the  correlator  is  controlled 
by  the  low  speed  logic,  in  particular,  the  select  pulse 
generator.  After  turn  on,  the  start  line  is  manually 
held  high  until  the  measurement  is  to  commence.  Thereafter, 
one  load  pulse  will  be  counted  before  the  first  select 
pulse  provides  an  interrupt  for  the  Nova,  and  the  buffer 
RAM  is  loaded  into  Nova  core.  Since  there  is  no  way  of 
determining  the  initial  medium  speed  counter  status,  this 
first  data  set  is  ignored.  The  following  full  cycle  will 
then  provide  the  first  good  data  set. 

9 . 2  Reference  Switching 

The  next  phase  is  to  view  a  reference  source, 
achieved  by  switching  the  frequency  of  the  radiometer 
local  oscillator.  This  will  provide  a  background  noise 
signal  equivalent  in  temperature  (see  Figure  1)  and  vari¬ 
ance  to  the  signal  frequency  band.  Thus,  viewing  this 
background  reference  for  an  equal  time  and  subtracting  this 
information  from  the  signal  information,  will  also  subtract 
out  any  signal  or  receiver  DC  bias.  The  result  will  be  a 
relative  measurement  with  only  signal  information.  This 
will  be  accomplished  by  the  receiver  providing  two  sets  of 
signals  shown  in  Figure  31.  Signal  A  is  applied  to  the 
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invert  control  of  Figure  25  and  controls  whether  the  data 
link  is  added  or  subtracted  from  the  running  sum  stored  in 
RAM.  Signal  B  is  applied  to  the  external  reset  line  on 
Figure  2S .  This  resets  the  medium  speed  counters  and  stalls 
the  low  speed  accumulator  until  the  radiometer  local 
oscillator  locks  up  on  a  new  frequency.  The  high  speed 
sections  have  been  running  constantly,  thus  when  Signal  B 
goes  low,  valid  data  has  already  been  accumulated  in  the 
high  speed  counters. 

Then  to  provide  information  pertaining  to  receiver 
gain  variations,  the  receiver  input  is  switched  from  the 
antenna  to  a  reference  source  via  a  waveguide  switch.  Fre¬ 
quency  switching  is  still  maintained  while  viewing  the 
reference,  thus,  the  reference  temperature  must  be  changed 
between  the  on  and  off  frequency  positions  in  order  not  to 
give  a  zero  result  as  a  reference.  The  reference  information 
is  to  be  stored  separately  from  the  signal  information, 
thus  the  low  accumulation  time  provides  a  good  division. 

This  would  make  2  x  96.5  seconds  the  minimum  switching  period. 
However,  since  a  50%  duty  cycle  between  the  antenna  and 
reference  source  may  not  be  desired,  a  status  bit  will  be 
set  in  order  to  notify  whether  the  incoming  data  is 
signal  or  reference.  This  circuitry  is  shown  in  Figure  32 
and  is  built  on  board  2  of  the  wire-wrap  series.  Therefore, 
when  the  radiometer  desires  a  reference  block,  it  sets  the 
reference  line  high  and  waits  for  a  select  pulse.  On  the 


Figure  31.  Reference  timing 
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Figure  32.  Status  generator 
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rising  edge  of  Che  select  pulse,  the  external  reset  line 
of  Figure  28  should  be  set  high.  The  wave  guide  switch  is 
then  switched  and  after  a  setting  time,  the  external  reset 
line  is  forced  low.  This  process  is  similar  to  the  nulling 
action  during  fi'equency  switching.  The  status  bit  then 
provides  the  status  for  the  next  data  set. 

The  next  step  is  to  read  sets  of  data  into  the  Nova. 
This  process  is  initiated  by  the  select  line  going  high  and 
interrupting  the  Nova.  Actual  storage  is  done  by  building 
a  memory  block  in  core  consisting  of  the  data  sets  accumu¬ 
lated  from  each  low  speed  count  period.  An  alternative  is 
to  accumulate  low  speed  data  sets  in  Nova  core,  that  is, 
one  data  set  would  be  recorded  and  additional  sets  would 
simply  be  added  to  the  existing  set.  Software  for  the 
former  case  is  given  in  Appendix  D.  Finally,  the  memory 
block  is  written  on  9  track  tape  for  final  reduction  by 
the  IBM  370  system. 

9 . 3  Quantization  Corrections 

The  accumulated  raw  data  are  now  divided  by  the 
number  of  samples  taken  to  give  the  final  "quantized"  auto¬ 
correlation  function.  The  number  of  low  speed  data  sets 
accumulated  is  the  first  word  of  the  memorv  block.  Each 
data  set  then  consists  of  16  high  speed  counts  times  2047 
medium  speed  counts.  For  a  96  second  correlation  time 
the  medium  speed  counter  is  counted  262,144  times,  thus, 

K  in  Equation  (6)  is  8 , 585 , 740 , 2S8  . 
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The  resultant  "quantized"  autocorrelation  function 
is  normalized  to  the  zero  delay  channel  and  corrected  for 
coarse  quantization  error.  For  the  zero  delay  channel, 
i=0  (Equation  13)  thus  r(iAx)  and  p (iAx) (normalized 
corrected  autocorrelation  function)  equals  1.  The  zero 
delay  channel  then  represents  the  total  power  received. 

For  other  channels,  i_>l ,  p(iAi)  is  much  less  than  1 
(Gaussian  input  uncorrelated)  thus  p(iAx)  is  closer  to  1. 
Equation  14  can  be  used  truncating  the  terms  in  Equations 
25,  27,  28,  and  29  to  provide  desired  accuracy. 

9 . 4  Weighting  and  Denormalization 

Next  the  normalized  autocorrelation  function 
p(iAx)  is  multiplied  by  a  weighting  function,  w(iAx) . 

As  previously  stated,  the  choice  of  w(iAx)  is  not  critical, 
and  three  alternatives  are  adapted  from  Weinreb  (1965)  in 
Table  2. 

A  normalized  power  spectrum,  p(iAf)  is  given  by 
the  Fourier  Transform  of  o(iAx).  The  system  gain  can  then 
be  determined  from  the  reference  measurement  and  incorpo¬ 
rated  into  the  signal  measurement  to  provide  o(iAf) 

(Equation  7) ,  the  actual  power  spectrum. 
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CHAPTER  X 

EXPERIMENT  DURATION 

10 . 1  System  Operating  Temperature 

The  length  of  time  during  which  the  experiment  is 
performed  is  based  on  the  minimum  detectable  signal  as 
given  in  Equation  (1).  In  this  case,  a=l  for  a  full 
power  receiver  (Tiuri,  1964),  and  BW  and  Top  need  to  be 
estimated . 

Top  for  a  double  side  band  receiver  where  a  signal 
appears  only  in  one  side  band  is  given  below  (adapted 
from  Okean  and  Lombardo,  1973): 

Top  =  2 (Ta+Te)  (40) 


Where : 

Ta  =  Temperature  seen  by  antenna 
Te  =  Equivalent  receiver  temperature 
Te  can  be  determined  by  considering  Figure  6.  Assuming 
the  waveguide  switch  loss  negligible,  Te  will  depend  pri¬ 
marily  on  the  mixer-preamplifier  noise  figure.  The  mixer- 
preamp  used  is  a  SpaceKom  model  F22.3V.  An  RF  to  IF  gain 
of  34  db  will  cause  noise  introduced  in  succeeding  stages 
to  be  negligible.  Thus,  the  mixer-preamp  broadband  double 
sideband  noise  figure  of  2.6  db  will  give  Te  from  the 
following  relationship  (Mumford  and  Scheibe,  1968): 

Te  =  (FDSB-1)290°K 
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Where : 

^DSB  =  double  sideband  noise  figure  with  signal 
present  in  both  bands  (Pure  Number") 

Te  is  then  equal  to  237.7  K.  Ta  can  be  approximated  from 
Figures  4  and  5  as  173.64  K  and  433b. 0  K  for  the  emission 
and  absorption  experiments,  respectively.  Thus,  Top  for 
emission  and  absorption  will  be  823  K  and  9147  K. 

10 . 2  Overall  Integration  Times 

With  respect  to  bandwidth,  the  system  operates  as 
the  equivalent  of  16  full  power  receivers  since  the  corre¬ 
lator  consists  of  16  channels  of  pure  integrators.  Thus, 

Tj.  will  be  considered  a  pure  integration  time,  and  RW 
will  be  approximated  as  20  MHz/ 16  or  1.25  MHz. 

Other  factors  influencing  r  j.  are  the  time  spent 
viewing  reference  sources  and  the  loss  encountered  by 
resetting  the  high  speed  accumulator.  Viewing  references 
will  reduce  by  1/4,  neglecting  the  relatively  short 
blanking  time  (Signal  B,  Figure  37").  Skipping  2  of  every 
18  samples  in  the  high  speed  accumulator  then  further  re¬ 
duces  the  actual  time  spent  correlating  to  8/9.  A  final 
factor  affecting  integration  time  is  the  quantization  error 
As  previously  stated,  for  the  oversampling  case,  is 
effectively  reduced  by  1/1.43. 

Bv  definition,  AT  .  of  Equation  (1)  is  derived  from 

mm  n 

the  idea  that  the  minimum  detectable  signal  is  equal  to  the 
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standard  deviation  of  the  received  power  (Tiuri,  1964). 

In  practice,  however,  the  minimum  detectable  signal  should 
be  made  1/4  of  the  received  power's  standard  deviation. 
This  new  value  of  AT  .  will  minimize  the  probabilitv  of 
random  noise  being  interpreted  as  signal  (Tiuri,  1964). 
However,  from  the  original  design  criteria,  signal  detec¬ 
tion  is  not  sufficient.  An  accurate  representation  of 
spectral  line  shape  was  originally  specified.  The  overall 
accuracy  will  now  be  defined  as  1%  of  the  signal  amplitude 
limited,  of  course,  to  the  16  equally  spaced  spectral 

points.  Thus,  AT  .  1 s  of  0.005  K  and  0.30  K,  for  the 
r  mm 

emission  and  absorption  experiments,  respectively  (from 
Figures  4  and  5) ,  will  provide  reasonable  minimum 
integration  times. 

Equation  (1)  now  becomes: 


AT. 


min 


Top 


[(1.25x106H2)(ti) (4) (9/8) (1.43)] 1/2 


(42) 


For  the  emission  experiment,  tj  becomes  approximately 
3400  seconds  and  for  absorption,  115  seconds. 


** _ 
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CHAPTER  XI 
CONCLUSIONS 


II. I  System  Choice 


In  this  study,  a  system  is  described  that  will  detect 
stratospheric  and  mesospheric  water  vapor.  The  method  of 
detection  involves  a  remote  measurement  of  the  molecular 
absorption  or  emission  spectrum.  The  resultant  spectral 
data  can  then  be  used  to  determine  water  vapor  mixing 
ratios  versus  altitude  by  inverting  the  radiative  transfer 
equation.  The  accuracy  of  this  type  of  measurement  is 
heavily  dependent  on  spectral  line  shape,  thus  much  of 
this  study  concerns  a  system  capable  of  accurate  spectral 
measurement . 

The  system  consists  of  a  22  GHz  receiver  feeding  a 
20  MHz  bandwidth  spectral  analyzer.  The  emphasis  then  is 
placed  on  the  design  of  a  suitable  spectral  analyzer  of 
which  two  types  were  considered.  These  were  a  discrete 
filter  bank  and  an  autocorrelator  in  both  analog  and 
digital  forms.  With  design  considerations  of  stability, 
sensitivity,  resolution,  complexity,  and  cost,  a  1  x  2 
bit  digital  autocorrelator  was  chosen  as  the  best  overall 
system. 


11.2  System  Design 

Sampling  at  the  Nyquist  rate  is  adequate  in  multibit 
systems.  However,  for  a  hard  quantization  scheme,  as  is  the 
case  here,  system  sensitivity  is  reduced  but  can  be  partly 


recovered  by  sampling  at  twice  the  Nyquist  rate  (80  Mhz) . 
Also,  a  further  increase  of  the  sampling  rate  to  100  MHz 
will  allow  adjustment  of  the  resultant  frequency  centers 
without  decreasing  sensitivity. 

The  autocorrelation  and  first  accumulation  are 
then  constructed  on  printed  circuit  boards  using  the  ECL 
10,000  logic  series.  A  medium  speed  accumulator  is  then 
also  constructed  on  printed  circuit  boards  using  both 
Schottky  and  standard  TTL  logic.  A  final  low  speed  accu¬ 
mulator  and  computer  interface  is  then  wire  wrapped  using 
MOS  3nd  TTL  logic.  All  further  data  reduction,  including 
a  discreet  Fourier  Transform  and  spectrum  denormalization, 
is  performed  in  software. 

11.3  Suggestions  for  Future  Work 

At  the  time  of  this  writing,  a  prototype  system 
consisting  of  16  autocorrelation  points  was  designed,  two 
points  of  which  were  built.  An  immediate  goal  then  is  to 
complete  the  16  points  and  use  the  system  for  a  ground 
based  measurement. 

To  complement  the  hardware  system,  a  software  pack¬ 
age  can  be  constructed  to  perform  various  tasks.  For 
instance,  different  types  of  Fourier  Transforms  (FFT, 
Chirp-Z)  can  be  programmed  to  allow  grouping  or  spreading 
of  spectral  points.  Also,  overall  teletype  control  and  rea 
time  displays  would  enhance  system  operation. 


Extended  goals  would  he  to  perform  an  intra-atmosphe 
measurement  from  a  balloon.  Such  a  system  would  require  a 
telemetry  link  to  ground  and  could  also  provide  real  time 
results.  A  third  alternative  is  a  "look  down"  experiment, 
that  is,  a  measurement  from  outside  the  earth's  atmosphere. 
This  type  of  measurement  might  be  performed  in  the  19S0's 
as  part  of  a  Space  Shuttle  experiment. 
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APPENDIX  A 
TRUNCATION  ERROR 

Quantifying  the  error  introduced  from  truncating 
counter  bits  can  be  handled  by  considering  the  effects  of 
variance  on  the  total  resultant  counts.  Abies,  et  al. 
(1974)  handled  a  similar  truncation  problem  by  adding  a 
"round  off"  variance  to  the  full  counter  variance.  Here 
then,  dividing  the  difference  between  the  full  counter 
plus  rounded  counter  variance  and  the  full  counter  variance 
by  the  full  counter  variance  will  provide  a  relative 
indication  of  the  change  in  overall  variance. 

2 

The  rounding  error  (variance)  is  given  by  * 

2 

V“/l 2  where  V  is  the  value  of  the  least  significant  bit 
used.  The  probability,  p,  of  counts  occurring  is  given  by 
the  numerator  of  Equation  (14) ,  used  in  conjunction  with 
Equations  (30)  through  (33)  .  The  variance  of  a  full  digit 
counter  is  then  Cp^  =  P(1-P)N,  where  N  is  the  total  number 
of  samples.  The  relative  variance  is: 


Using  Vo-0.95,  N-4,  and  V-23 : 


(A3) 


2  _  23.02 

'  REL  "  nr 

Obviously  then,  for  loosely  correlated  signals,  as 

is  the  case  here,  N  will  need  to  be  extremely  large  in  order 

to  provide  an  insignificant  change  in  system  sensitivity  due 

to  bit  truncation.  Thus,  if  the  incoming  signal  is 

assumed  very  loosely  correlated,  say  paO.OOOl,  then  an  N 

8 

on  the  order  of  10  is  necessary  to  reduce  the  change  in 

g 

variance  to  0.17o.  However,  at  a  100  MHz  sampling  r3te,  10 
samples  occur  every  second.  Therefore,  the  truncation 
error  can  be  ignored  for  correlation  time  on  the  order 
of  seconds  or  longer. 
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APPENDIX  B 
TIMING  DIAGRAMS 

Following  are  expected  waveforms  from  throughout 
the  system.  Note  these  waveforms  are  of  two  types,  one 
describing  actual  signal  levels  and  the  other  describing 
whether  data  are  valid  or  not.  The  "data  valid"  waveforms 
will  be  labeled  in  parenthesis  and  will  indicate  data 
valid  as  high,  invalid  as  low.  The  signal  level  waveforms 
are  labeled  without  parenthesis  and  represent  a  logical 
one  when  high,  a  zero  when  low. 


1  • 


4k  i  INI  I 


Figinre  A2.  High  spec-d  data  propagation 


AITKNDIX  C 


WIRE  WRAP  HOARDS 

Shown  hero  are  the  physical  layouts  o t  the  two 
wire  wrap  boards .  hoard  1  consists  of  t.ho  low  speed 
accumulator  as  shown  in  Figure  2o  (top)  whereas  Board  2 
consists  predominantly  of  computer  interfacing  circuitry 
(bottom  of  Figure  26)  and  bookkeeping  (Figures  .10  and  .11). 

At  the  time  of  this  printing,  no  radiometer  inter¬ 
facing  circuitry  was  constructed,  although  ample  space 
was  available  on  either  wire  wrap  board.  One  such  inter¬ 
face  was  presented  in  Figure  33. 

Also  included  here  are  the  pin  connections  for  the 
two  wire  wrap  boards  and  the  bus  connections  to  the  Data 
General  Nova  1200  minicomputer. 


jlA 
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NOVA  RACK  FRONT  VIEW 
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\b  input/ 
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CONNECTOR  FRONT  VIEW 
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Figure  A7.  Nova  data  interface  (after  Sweitzer,  1978) 
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Wire  wrap  Board  L,  Connector  A 
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Wire  wrap  Board  1,  Connector  B 
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APPENDIX  D 
NOVA  SOFTWARE 

Presented  here  is  the  software  routine  enabling 
data  to  be  stored  in  the  Data  General  Nova  1200  minicomputer. 
Data  are  stored  in  the  form  of  32  bit  words,  forming  a 
memory  block,  beginning  at  location  200g.  The  higher  order 

16  bits  of  a  word  are  stored  first  (lower  memory  location) . 
The  first  word  of  the  memory  block  provides  the  number  of 

17  word  cycles  recorded  in  the  block.  Next,  the  first  word 
of  a  cycle  is  recorded,  providing  the  status  of  that  cycle. 

A  1073741824  (bit  2^  set)  indicates  signal  data  and  a 
3221225472  (bits  2^  and  2"^  set)  indicates  reference 
data.  Finally,  16  data  words  are  recorded,  beginning  with 
zero  lag  point.  The  17  word  cycles  are  then  repeated  until 
termination  of  the  experiment. 

The  program  operates  as  an  interrupt  service  routine, 
that  is,  the  program  sits  in  a  loop  in  between  the  17  word 
cycle  recordings.  A  select  pulse  then  interrupts  the  loop 
and  a  cycle  of  data  is  recorded.  Upon  reception  of  a  zero 
status  word,  the  program  is  terminated. 

For  transfer  to  the  IBM  370  system,  the  Nova  can 
record  the  data  block  directly  on  9  track  tape,  location 
2 00 q  first.  Then,  using  Fortran  IV,  the  IBM  system  can 
directly  form  10  bit  (14)  integers. 
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Interrupt  Service 
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up  number  of  words  to  read  on  tape 
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